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1.  INTRODUCTION 


The  probe  reported  here  was  designed  to  be  used  to  deter¬ 
mine  the  radial  distribution  of  the  time-averaged  flow  in  a  small, 
single-stage  axial  compressor  (Figure  1) .  While  the  probe  has 
more  general  application,  the  chosen  geometry  is  the  result  of 
its  first  intended  application. 

The  probe  was  designed  to  measure  air  flows  in  the  range  of 
Mach  number  from  0.3  to  0.7  with  pitch  angles  expected  in  the 
range  0  to  15  degrees  and  stagnation  temperatures  from  atmospheric 
to  about  150°F  (65°C) .  Static  pressure  would  range  from  0.6  to 
about  1.2  atmospheres. 

The  prototype  of  the  present  probe  was  built  by  Dodge  and 
was  reported  in  Ref  1.  The  first  application  of  the  earlier 
design  to  measure  rotor  losses  and  deviation  angles  was  reported 
in  Ref  2  and  Ref  3.  In  Ref  3,  errors  in  the  earlier  work  were 
corrected  and  the  method  used  to  analytically  represent  the  cali¬ 
bration  of  the  probe  was  set  out  more  clearly. 

The  design  of  the  present  probe  departs  only  in  detail  from 
that  of  the  prototype,  for  reasons  which  are  given  in  Section  2. 
However,  the  method  now  used  to  represent  the  calibration  of  the 
probe  is  quite  different  and  considered  to  be  much  improved  over 
that  reported  in  earlier  references. 

In  the  present  report,  the  probe  design  is  described  in 
Section  2  and  the  analytical  and  experimental  procedures  used  in 
its  calibration  are  given  in  Section  3.  Steady  flow  tests  to 
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verify  the  calibration  and  compressor  measurements  are  reported 
in  Section  4  and  conclusions  are  given  in  Section  5.  The 
analytical  basis  for  the  pneumatic  calibration  representation 
is  given  in  Appendix  A.  The  computer  programs,  procedures  and 
data  files  used  in  the  calibration  and  verification  are  given 
in  Appendix  S. 
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2.  PROBE  DESCRIPTION 


2.1.  General  Arrangement 

The  probe  is  designed  to  measure  four  pressures  and  one 
temperature.  Figure  2  shows  a  drawing  of  the  probe  and  a  view 
is  shown  in  Figure  3.  Due  to  the  symmetry  of  the  probe  about 
one  plane  through  its  (radial)  axis,  a  balancing  of  the  pres¬ 
sures  at  sensors  P2  and  P3  in  Figure  3  by  rotating  the  probe 
around  the  radial  axis  allows  the  yaw  angle  to  be  read  against 
a  vernier  scale  or  recorded  using  a  potentiometer  read-out  of 
the  angular  orientation.  The  pitch  angle  as  well  as  the  mag¬ 
nitude  of  the  dimensionless  velocity  (or  Mach  number)  at  the 
probe  tip  are  calculated  using  the  pressure  readings  of  the 
four  sensors.  The  relationship  of  the  pressures  to  Mach  number 
and  pitch  angle  must  first  be  established  by  calibration.  The 
stagnation  temperature  is  needed  in  order  to  calculate  the  mag¬ 
nitude  of  the  velocity  from  the  dimensionless  velocity  or  Mach 
number.  The  exposed  fine-wire  thermocouple  sensor  supported  in 
the  stagnation  region  formed  by  a  glassy  aluminum  oxide  insula¬ 
tor.-  gives  an  output  which  corresponds  fairly  closely  to  stag¬ 
nation  temperature.  The  departure  from  stagnation  temperature, 
the  recovery  factor ,  must  also  be  established  by  calibration. 

2.2.  Design  Features 

The  probe  was  designed  to  have  separate  sensors  to  avoid 
the  so-called  "stem-effect"  associated  with  cylindrical  probes 
having  surface  sensors  (Reference  4) .  The  restricted  axial  gap 
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in  the  compressor  (0.63")  did  not  allow  the  use  of  a  "goose¬ 
neck"  conical  probe.  As  was  reported  in  Reference  (2) ,  the 
stem-effect  was  found  to  be  negligible  for  the  first  probe 
geometry  to  within  0.25"  of  the  wall.  The  probe  was  designed 
to  be  translated  through  a  hole  of  0.25  inches  diameter.  The 
sensors  PI,  P2  and  P3  were  placed  at  the  same  radial  location 
and  therefore  similarly  average  the  flow  from  the  upstream 
rotor  blade  wakes.  These  sensors  largely  dictate  the  measure¬ 
ment  of  Mach  number.  The  pressures  PI  and  P4  largely  dictate 
the  measurement  of  pitch.  The  sensor  P4  was  unavoidably  dis¬ 
placed  radially  from  the  other  three.  However,  since  the  pitch 
angle  variation  to  be  measured  was  rather  small,  and  since  the 
radial  gradients  in  flow  properties  were  small  outside  the  wall 
boundary- layers ,  this  necessity  was  easily  accepted.  It  is 
also  noted  that  when  the  probe  is  adjusted  in  yaw,  the  sensor 
PI  will  indicate  stagnation  pressure  very  closely  since  the 
pitch  angle  is  not  expected  to  exceed  15°. 

The  radial  separation  of  the  temperature  sensor  was  also 
easily  accepted  since  the  probe  was  to  be  translated  in  radial 
surveys  and  the  physical  displacement  could  be  accounted  for 
in  the  data  reduction.  The  arrangement  and  the  geometry  of 
the  temperature  sensor  were  chosen  so  that  a  high  recovery 
factor  and  sensitivity  would  be  obtained  and  yet  the  element 
could  be  repaired  easily  in  the  event  of  a  failure.  The  shape 
ground  into  the  glassy  insulator  and  the  method  of  stringing 
the  thermocouple  wire  were  changed  from  those  reported  in  Ref¬ 
erences  1  and  3.  In  the  new  arrangement,  a  long  length  of 
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thermocouple  wire  is  exposed  on  either  side  of  the  junction 
within  the  near-stagnation  flow  ahead  of  the  insulator.  This 
results  in  a  higher  recovery  factor  at  all  Mach  numbers  and  a 
more  easily  controlled  and  repeatable  geometry. 

The  tubes  for  sensors  P2,  P3  and  P4  (Figure  3)  were 
oriented  at  55°  to  the  sensor  PI  based  on  the  known  behavior 
of  a  straight  cylindrical  tube  when  pitched  or  yawed  to  the 
flow  (Reference  1) .  An  angle  of  about  55°  was  required  in  order 
for  the  pressure  at  the  sensor  to  decrease  from  stagnation  to 
static  pressure.  In  the  earlier  design  (Reference  1)  the  centers 
of  the  sensors  were  arranged  to  be  on  the  surface  of  a  small 
sphere  which,  in  effect,  determined  the  spatial  resolution  of 
the  probe.  This  is  shown  in  Figure  4.  Also  in  Figure  4  is 
shown  the  tip  design  of  the  present  probe.  The  sensors  were 
initially  positioned  such  that  the  centers  of  the  tubes  were 
aligned  in  a  plane,  and  were  closer  together  than  in  the  earlier 
design.  Initial  tests  of  the  revised  geometry  resulted  in  the 
measurement  of  small  differences  between  the  pressure  sensed  at 
P2,  P3  and  P4,  and  that  at  Pi.  It  was  conjectured  that  the 
surrounding  sensors  were  too  nearly  aligned  with  the  flow  in¬ 
duced  over  the  central  sensor  and  were  therefore  measuring  close 
to  stagnation  pressure.  The  tips  of  the  surrounding  sensors 
were  modified  as  shown  in  Figure  4  and  pressure  differences 
near  to  dynamic  pressure  were  subsequently  obtained. 
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3 .  CALIBRATION 


3.1.  Apparatus  and  Procedure 

The  probe  was  calibrated  in  a  4.25"  diameter  free jet  ex¬ 
hausting  to  atmosphere.  The  apparatus  shown  in  Figure  5  can  be 
used  to  establish  air  flows  to  a  Mach  number  of  unity.  Pitch 
angle  can  be  set  within  a  range  of  +45  to  -45  degrees  while 
the  yaw  angle  can  be  set  to  any  angle  desired.  In  effecting 
angle  changes  the  probe  tip  is  rotated  about  its  tip  which 
remains  at  the  same  point  on  the  center  line  of  the  jet. 

The  geometry  for  the  tests  and  the  instrumentation  are 
shown  in  Figure  6.  The  speed  of  the  flow  was  adjusted  by  moni¬ 
toring  the  impact  pressure  using  a  water  manometer  board.  The 
static  pressure  was  taken  to  be  the  atmospheric  pressure  while 
the  total  pressure  and  temperature  were  measured  in  the  air 
supply  pipe  immediately  upstream  of  the  nozzle.  No  difference 
in  total  pressure  was  detected  in  the  flow  between  the  sensor 
upstream  of  the  nozzle  to  the  core  of  the  jet  and  hence  the 
same  total  conditions  (pressure  and  temperature)  were  assumed 
for  the  pipe  and  the  jet  core.  Since  the  temperature  of  the 
supplied  air  fluctuated  measurably  due  to  ambient  (coastal) 
variations,  the  difference  between  the  total  temperature  in  the 
pipe  and  the  temperature  detected  by  the  combination  probe  was 
recorded  differentially.  A  thermocouple  probe  having  a  fine- 
wire  sensor  similar  to  that  of  the  combination  probe  was  used 
at  the  upstream  station.  The  absolute  value  of  the  temperature 


in  the  supply  air  was  measured  using  the  same  upstream  thermo¬ 
couple  element  referenced  to  ice  point. 

All  pressures  were  converted  to  voltages  using  a  Scani- 
valve  and  single  calibrated  transducer  for  which  the  reference 
pressure  was  atmospheric.  The  thermocouple  voltage  differences 
were  recorded  without  amplif ication.  Yaw  angles  were  adjusted 
to  a  vernier  scale  and  converted  to  voltage  using  a  potentio¬ 
meter.  The  pitch  angle  was  set  to  a  scale  on  the  adjustable 
mount  and  recorded  manually. 

Pressure  tubes  P2  and  P3  were  connected  to  a  U-tube  mano¬ 
meter.  The  probe  was  aligned  (in  yaw)  with  the  flow  by  balancing 
the  two  pressures.  The  pneumatically  averaged  value  of  the 
two  pressures  -  designated  P23  -  was  recorded  and  used  in  repre¬ 
senting  the  calibration  of  the  probe.  All  data  were  recorded 
and  analyzed  using  the  data  acquisition  system  shown  in  Figure 
7.  The  programs  and  organization  of  data  files  are  given  in 
Appendix  e.  Before  beginning  calibration  tests,  since  the 
sensors  for  two  pressures  and  the  sensor  for  temperature  were 
physically  separated  radially,  the  velocity  distribution  across 
the  jet  was  examined.  Figure  8  shows  the  velocity  distribution 
obtained  by  probe  measurement.  It  can  be  seen  that  for  a  core 
of  approximately  3  inches  in  diameter  the  flow  was  uniform. 

The  procedure  was  as  follows: 

The  required  Mach  number  was  established.  Having  allowed 
some  5  to  10  minutes  for  the  flow  to  stabilize,  the  yaw  angle 
of  the  probe  at  zero  degrees  of  pitch  was  checked.  With  the 
probe  aligned  with  the  flow,  P2  and  P3  should  be  the  same  and 
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the  angle  read  on  the  vernier  scale  should  be  zero.  Any  depar¬ 
ture  was  corrected  by  adjusting  the  probe  in  the  probe  holder. 
The  first  pitch  angle  was  set  at  -6°  and  a  scan  of  the  data  was 
taken.  The  pitch  angle  was  changed  in  increments  of  2°  up  to 
an  angle  of  +18°,  taking  a  data  scan  at  each  setting.  A  sample 
of  the  raw  data  output  is  shown  in  Table  I. 

After  one  survey  over  pitch  angle,  the  probe  was  reset  to 
zero  pitch  and  zero  yaw,  the  Mach  number  of  the  flow  was  changed 
and  the  procedure  repeated.  The  range  of  Mach  number  from  0.3 
to  0.7  was  covered  in  increments  of  0.1  or  less. 


3.2.  Dimensionless  Velocity  -  Pneumatic  Calibration 

The  characteristics  of  the  probe  must  be  represented  such 
that  pressures  and  temperatures  measured  in  a  flow  can  be  related 
to  the  pitch  and  yaw  angles  and  the  velocity  magnitude.  How¬ 
ever,  since  the  probe  is  always  rotated  to  balance  the  pressures 
P2  and  P3,  the  yaw  angle  can  always  be  read  directly.  It  is 
left  therefore  to  establish  the  relationship  for  pitch  angle 
and  velocity  magnitude  in  terms  of  the  probe  pressures  and  tem¬ 
perature  for  the  special  case  of  zero  yaw. 

Appendix  A  discusses  the  reduction  of  the  four  measured 
pressures  to  basically  two  coefficients: 


-  _  PI  -  P23 

0  *  - pi - 


(1) 


and 


_  PI  - 

PI  -  P23  ' 


(2) 
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where 


5  =  y  *  S  (3) 

might  be  looked  at  as  an  alternate  choice.  These  coefficients 
can  be  established  using  the  pressure  readings  from  the  probe 
for  any  flow  condition.  The  calibration  gives  a  total  of  104 
different  values  for  0 ,  and  corresponding  values  of  y  and  <5 
for  104  separate  combinations  of  Mach  number  and  pitch  angle. 

The  reduced  data  are  given  in  Table  II. 

It  can  be  argued  (Appendix  A) ,  and  it  can  be  seen  in  the 
results  that  0  represents  largely  a  measurement  of  the  Mach 
number  and  y  provides  mainly  a  measurement  of  the  pitch  angle. 
However,  if  explicit  relationships  for  the  reduced  velocity,  X 
(defined  by  Equation  A(3)  using  Equation  A(l)),  and  the  pitch 
angle,  <p ,  in  terms  of  0  and  y  can  be  obtained  by  representing 
the  calibration  data  using  monotonic  mathematical  functions,  X 
or  $  can  then  be  determined  uniquely  for  any  given  values  of  0 
and  y.  The  method  involves  using  mathematical  approximations 
of  surfaces  which  represent  the  calibration  data.  In  the  way 
that  data  points  which  depend  on  just  one  variable  can  be  ap¬ 
proximated  by  a  polynomial,  X  and  $  are  approximated  as  being 
polynomial  functions  of  both  3  and  y.  This  leads  to  mathematical 
expressions  for  the  calibration  surfaces  which  are  of  the  form: 

x  «  l  {  l  }  *  yU-1) 

i-1  j-1  13 

and 


IS 


(4) 


4»  -  I  {  ?Di.0(j“1)}  •  y(l"1)  (5) 

i»l  j=l 

where  and  are  coefficients  which  must  be  derived  from 
the  data. 

Using  the  computer  programs  given  in  Reference  5,  the 
surfaces  represented  by  Equation  (4)  and  Equation  (5)  were  de¬ 
rived  from  the  calibration  data.  Also  derived  and  examined 
were  the  alternate  surfaces  obtained  by  electing  to  use  the 
coefficients  6  and  y  instead  of  0  and  y,  related  through  Equa¬ 
tion  (3) .  The  equations  for  the  alternate  surfaces  are  the 
following: 

X  *  l  {  Ic..a  }  •  Y(l_1)  (4a) 

i=l  j-113 

=*  l  {  ?D..<$(j_1)}  •  y11'11  (5a) 

i»l  j^l1-1 

[It  is  noted  that  the  representation  achieved  using  Equa¬ 
tion  (4)  and  Equation  (5)  or  Equation  (4a)  and  Equation  (5a)  is 
similar  to  the  earlier  method  of  polynomial  approximation 
described  in  Reference  3.  However,  in  the  method  of  Reference 
3,  the  expression  for  X  was  implicit  and  an  iterative  technique 
was  required  to  obtain  X  and  <J>  from  measured  values  of  8  and  y 
or  5  and  y.  The  derivation  of  the  polynomial  coefficients 
using  matrix  operations  makes  possible  the  use  of  explicit 
relations  for  both  X  and  $ . ] 

The  choice  of  using  Equation  (4)  and  Equation  (5)  or  Equa¬ 
tion  (4a)  and  Equation  (5a)  was  made  after  an  examination  of  the 
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relative  accuracy  of  the  two  methods.  Using  the  calibration 
data,  for  each  given  pair  of  values  for  <S  and  y,  X  and  <t>  were 
calculated  using  the  coefficients  determined  in  the  two  sets 
of  equations.  The  calculated  X  was  compared  in  each  case  to 
the  actual  value  set  in  the  calibration  test  and  an  error  was 
defined  as  follows: 


where  X.  is  the  actual  value  in  the  calibration  test  and  X„  is 
0  c 

the  value  calculated  using  the  surface  approximation  equation. 

The  error  in  X  so  defined  is  then  a  percentage  deviation  from 
the  actual  value.  For  the  pitch  angle,  the  error  as  a  percentage 
is  not  meaningful  (for  example,  when  3  0)  .  The  error  was 
defined  as 

e*  *  *0  “  (7) 
where  ^  is  the  actual  value  set  in  the  calibration  test  and  4>c 
is  the  value  of  the  pitch  angle  (in  degrees)  calculated  using 
the  derived  surface  approximation  equation.  From  an  examination 
of  the  relative  errors,  the  selection  of  using  Equation  (4a) 
and  Equation  (5a),  rather  that  Equation  (4)  and  Equation  (5), 
was  made. 

The  "surfaces"  obtained  by  drawing  straight  lines  between 
data  points  on  a  3D  plot  of  the  test  data  are  shown  in  Figures  9 
and  10  together  with  the  "surfaces"  obtained  by  joining  the  points 
calculated  using  Equation  (4a)  and  Equation  (5a) .  The  magnitude 
of  the  errors  shown  graphically  in  Figure  9  are  given  in  Table 


Ill,  together  with  the  values  of  the  coefficients  used  in  the 
approximation  equations . 


In  the  method  given  in  Reference  5  the  order  of  the  poly¬ 
nomial  approximation  can  be  changed  between  1  and  6  for  the 
dependence  on  S  and  between  1  and  6  for  the  dependence  on  y , 
independently.  The  plots  of  Figure  9  and  the  corresponding 
errors  in  Table  III  are  for  the  selection  which  gave  the  best 
results  in  that  the  averaged  error  was  least  in  the  range  of 
velocity  and  pitch  angle  anticipated  in  the  probe  application. 

As  indicated  in  Table  III,  the  range  expected  in  Mach  number  is 
from  .3  to  .7  and  in  pitch  angle,  from  0.0°  to  12.0°.  The 
maximum  error  within  this  range  was  +1.109%  in  X  and  0.417°  in 
pitch  angle,  while  the  average  error  was  -0.061%  in  X  and  -0.016° 
in  pitch  angle. 

3.3.  Temperature  Recovery  Calibration 

Since  the  probe  was  to  be  used  to  measure  losses,  requiring 
an  evaluation  of  velocity  magnitude  rather  than  simply  the  Mach 
number  or  dimensionless  velocity,  it  was  also  required  to  measure 
the  local  flow  total  temperature. 

Even  in  a  flow  in  which  the  angle  of  the  velocity  vector 
is  uniform  and  known  it  is  difficult  to  design  a  probe  which 
will  measure  true  stagnation  temperature  to  an  acceptable  ac¬ 
curacy  when  the  Mach  number  is  in  the  higher  subsonic  range. 

In  the  present  design,  advantage  was  taken  of  the  fine  wire 
sensor  being  small  and  in  poor  thermal  contact  with  a  glassy 
insulator.  While  the  "temperature  recovery"  of  the  probe  was 


expected  to  be  comparatively  high,  a  calibration  to  establish 
the  value  of  the  temperature  recovery  factor  was  required  never¬ 
theless.  In  the  calibration  test  a  voltage  was  recorded  as  the 
result  of  the  temperature  difference  between  the  combination 
probe  sensor  and  the  total  temperature  sensor  in  the  pipe  up¬ 
stream  of  the  jet  nozzle.  (The  relation  between  the  output  in 
millivolts  and  the  temperature  in  degrees  Fahrenheit  for  a  single 
sensor  referenced  to  ice  point  is  given  in  Appendix  C.)  The 
output  of  the  upstream  probe  was  also  recorded  with  reference 
to  ice  point.  As  the  upstream  probe  was  always  in  a  low  Mach 
number  flow  (less  than  0.2),  it  was  assumed  that  its  output 
corresponded  to  the  total  temperature  of  the  flow.  The  voltage 
difference  between  the  upstream  probe  and  the  combination  probe 
was  then  a  measure  of  the  departure  of  the  combination  probe 
temperature  from  the  stagnation  temperature. 

The  "temperature  recovery  factor",  R,  is  defined  as 


where  Tp  is  the  temperature  indicated  by  the  probe,  Tg  is  the 
static  temperature,  and  Tfc  is  the  total  temperature  of  the  flow. 

Since  Equation  (8)  contains  the  static  as  well  as  the 
total  temperature,  one  of  the  two  can  be  expressed  in  terms  of 
the  other  and  the  Mach  number  or  dimensionless  velocity.  Using 
Equation  A (4)  and  rearranging.  Equation  (8)  becomes 


- 1  "  '2  ' 

l  -  +  x‘r 


(9) 


where  X  is  the  velocity  made  dimensionless  by  dividing  by  the 
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"limiting”  or  "stagnation"  velocity  defined  by  Equation  A(l) . 

Since  X  is  established  by  applying  the  calibration  in  Section 
3.2  to  the  pneumatic  measurements,  it  is  only  necessary  to  know 
R  in  order  to  determine  from  the  temperature  Tp  indicated  by 
the  probe.  From  the  physical  arrangement  of  the  probe  it  is 
evident  that  the  thermocouple  junction  is  differently  exposed 
to  the  flow  at  different  yaw  and  pitch  angles,  and  the  Mach 
number  is  known  to  have  a  major  influence  on  the  temperature 
recovery  factor.  However,  since  the  yaw  angle  is  always  ad¬ 
justed  to  zero  and  the  pitch  angle  and  Mach  number  depend  uniquely 
on  the  measured  values  of  6  and  y,  the  recovery  factor  can  also 
be  approximated  as  a  surface  of  values  depending  on  6  and  y; 
thus 


M 


R 


[  {  I  Ei  S'5'1') 

i=l  j=l  13 


(i-1) 


(10) 


The  calibration  surface  obtained  from  the  calibration  test 
data  is  shown  in  Figure  11.  The  coefficients  evaluated  from  the 
data  and  the  resulting  percentage  errors  in  the  approximation 
are  given  in  Table  IV.  It  can  be  seen  that  errors  of  less  than 
0.8%  occur  within  the  range  of  measurements. 
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4.  APPLICATION 


4.1.  Verification  of  the  Calibration 

Before  using  the  probe  for  compressor  measurements,  the 
probe  was  mounted  on  the  free  jet  again  to  obtain  data  with 
which  to  verify  the  calibration.  The  probe  was  aligned  with 
the  flow  but  Mach  number  and  pitch  angle  were  set  remotely  and 
unknown  to  the  operator.  A  total  of  12  separate  conditions 
with  various  combinations  of  pitch  angle  and  Mach  number  were 
set.  The  results  are  shown  in  Table  V.  The  largest  error  was 

found  to  be  0.55°  in  pitch  and  0.7%  in  Mach  number. 

4.2.  Compressor  Measurements 

The  probe  can  be  used  in  any  flow  for  which  conditions 
are  within  the  ranges  of  pitch  angle  and  Mach  number  covered  in 

the  calibration.  Its  intended  application,  however,  was  to  de¬ 

termine  the  flow  between  blade  rows  in  the  compressor  shown  in 
Figure  1.  In  the  compressor  annulus,  large  radial  gradients 
might  be  expected  to  occur  (under  some  operating  conditions)  in 
pitch  angle,  yaw  angle  and  Mach  number,  so  that  account  had  to 
be  taken  of  the  physical  separation  (0.08  inches)  of  the  probe 
pressure  and  temperature  sensors  in  the  radial  direction.  This 
could  be  done  easily  since  the  probe  was  to  be  used  in  radial 
surveys  and  spatial  interpolation  could  be  used  in  reducing  the 
data. 

The  procedure  adopted  was  to  first  interpolate  the  raw 
data  obtained  for  the  radial  distribution  of  the  temperature 
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rise.  The  data,  in  the  form  of  the  differential  voltage  between 
the  probe  thermocouple  sensor  and  a  similar  reference  thermo¬ 
couple  sensor  in  the  compressor  inlet,  was  interpolated  using 
an  overlapping  quadratic  technique  to  obtain  the  differential 
voltage  at  the  locations  at  which  pneumatic  data  were  taken. 

The  temperature  recovery  factor  was  evaluated  at  each  location 
using  Equation  10  using  the  locally  measured  values  of  <S  and  y. 

The  probe  temperature  Tp  was  obtained  using  the  measured  value 
of  the  reference  temperature,  the  interpolated  differential 
voltage  and  the  thermocouple  calibration  given  in  Appendix  C. 

The  stagnation  temperature  was  obtained  using  Equation  9. 

An  example  of  the  radial  distribution  of  the  temperature 
rise  across  the  rotor  measured  in  this  way  is  shown  in  Figure 
12.  The  corresponding  data  are  given  in  Table  vi. 

The  uncertainties  in  the  measurements  made  in  the  com¬ 
pressor  flow  field  are  similar  to  the  uncertainties  present  in 
the  probe  calibration  measurements,  since  the  instrumentation 
used  was  the  same.  However,  the  calibration  was  carried  out  in 
a  uniform  steady  flow  whereas  the  flow  field  at  the  rotor  exit 
is  periodic,  therefore  unsteady.  The  possible  error  resulting 
from  unsteady  effects  is  discussed  in  Appendix  D.  For  the  data 
presented  here  the  possible  error  due  to  unsteady  effects  is 
negligible.  A  further  evaluation  of  unsteady  effects  must  be 
made  for  measurements  made  at  higher  rotational  speeds  and  flow 
Mach  numbers. 
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5.  CONCLUSIONS 


The  combination  probe  and  the  method  used  to  represent 
its  calibration  which  are  described  in  this  report  represent 
significant  improvements  over  the  prototype  probe  design  and 
calibration  method  reported  earlier.  The  pneumatic  characteris¬ 
tics  of  the  probe  were  improved  as  a  result  of  modifications  to 
the  tip  geometry.  This  allowed  the  pneumatic  characteristics 
to  be  well  ^represented  by  a  polynomial  surface  approximation. 

The  polynomial  surface  approximation  technique  used  to 
represent  both  the  pneumatic  and  temperature  characteristics  of 
the  probe  provides  a  simpler  and  more  accurate  technique  than 
was  previously  available.  Most  importantly,  however,  it  allows 
the  calibration  to  be  applied  directly,  allowing  Mach  number 
and  pitch  angle  to  be  calculated  from  pressure  measurements 
without  the  need  for  iteration. 

The  improved  probe  and  procedures  were  well  proven  in 
tests  carried  out  in  the  compressor. 
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102,483950 
102,348510 
102.379780 
192.386720 
112.435331 


T  ref  pipe 
93.913498 
99.083435 
99,270951 
99.492749 
99.546265 
99.745392 
99.967941 

144.199744 

'MM 

109.693869 
100.386980 
141 ,315399 


T  ref  pipe 
193.929814 
143.327360 
1*3.941244 
193,973364 
1*3.97626* 
144.99013* 
194.9445*4 
194.963844 
104.125151 
104.08720* 
104.13974* 
194.116410 
1*4,119311 


T  rtf  pipe 
104.256S04 
104,346129 

Mil 

i §4 , 32655# 
104,276934 
194.358661 
104.390769 
104.36742* 
104.288699 
194.33824* 
194.306121 
114.349151 


Tc»nbor»t>» 

95.592697 

95.968597 

94.394931 

94.296875 

95.016922 

94.9944*9 

94.623123 

95. 373281 
95.4698*9 
96.649658 
95.636749 
96.226959 
94.573166 


Tcw*pr»*« 

103.564589 
1*3.666789 
1*3.856570 
103.75436* 
103.799419 
t*3. 914959 
1*3.827360 
1*3.891629 
113.34196* 
143.914380 
104.0*546* 
193.77774* 
113.82736* 


Tc«nbor»l>e 

104.148519 

194.058119 

104.113469 

194.168959 

194.23314* 

114.148511 

1*4.095961 

104,168950 

104,12805* 

194.12S1SI 

194.020070 

104,093851 

193.891621 


Her.  Far 
.308953 
.321919 
.721063 
.705939 
.741226 
.o75W 
.694283 

.73*528 

:#8S 

.708847 

.734920 

.612831 


iec.Fnc 

.985369 

.993452 

.996525 

.991025 

.992332 

.992822 

.991254 

.992914 

.988412 

.992344 

.99452* 

.986105 

.98M69 


Her. Far 
,996016 
.99*835 
.995563 
.994381 
.996554 
.995267 
.990331 
,99182'’ 
.991149 
,993990 
.988262 
.992184 
.984545 


Tabic  II.  Complete  Got  of  Reduced  Calibration  Data 
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Oe*o 

.(53507 

.151745 

.(50401 

,149744 

,151291 

,150719 

.151715 

.152329 

:!SS) 

.153537 
, 153517 
.  153123 


Beta 
.173106 
.171937 
■  168925 
.167226 
. (68445 
.169123 
.171305 
.170497 
.170469 
.171751 
.172328 
.172361 
.172235 


Beta 

.194665 

.191932 

.19(324 

.199419 

.19(124 

.191159 

.190684 

.190860 

.192160 

.192600 

.193409 

.193438 

.194993 


Beta 

.217339 

.213911 

.210524 

.208658 

.212254 

.210802 

.210208 

.210208 

.210644 

.212185 

.212827 

.214413 

.214390 


Gown 

Delta 

i  tel 

mchna. 

1.595331 

.244894 

.240619 

.  554.528 

1.538058 

.233393 

.248285 

.553511 

( , 473545 

.221622 

.240397 

,55.5783 

1.394422 

.208807 

.248619 

.240508 

.554329 

1.300662 

.195478 

.554055 

1.208719 

.182177 

.240174 

.553238 

1.119423 

.169833 

.240731 

,554600 

1.043081 

.158892 

.240397 

.553783 

.959530 

.984416 

:186! 

■m 

,55.551 1 
.552692 

.822309 

.126255 

.240397 

.553783 

.767834 

.117868 

240285 

.5535(0 

.714844 

.109459 

.23983? 

.SS2419 

Gama 

Delta 

X  tel 

Hnchna. 

1.609(31 

.278389 

.258473 

.598294 

1.550336 

.266560 

.258869 

,599277 

(.484055 

.2S0694 

.2S8770 

.599031 

1.411968 

.236117 

.257677 

,596320 

1.3(5428 

.221577 

.2S8572 

,598540 

1.217540 

.205792 

.258473 

,598294 

1.126554 

.191857 

.2S7876 

.596814 

1.848533 

.178772 

.257677 

,596320 

.967232 

.164883 

.2S7876 

.596814 

.890134 

.823464 

.1S2881 

.141906 

m 

,598539 

.599031 

.773184 

.133267 

.258671 

.598786 

.717002 

,123493 

.257976 

.597061 

Goiwf 

Delta 

X  «al 

tactma. 

1.641706 

.279619 

. 651225 

1.583654 

.279533 

.651006 

1.506256 

.288183 

.279793 

.65(663 

1.427875 

.270466 

.279706 

.651444 

1.329150 

.254032 

.279880 

.65(883 

1.232625 

1.146035 

.235627 

.218531 

,278923 

,279359 

.649467 

.650566 

1.858027 

.201935 

.278662 

.648806 

.973129 

.186996 

.278574 

.648585 

.897510 

.172860 

.279098 

.649907 

.835878 

.161666 

.279793 

.65(663 

.779370 

.150760 

.278662 

.648806 

.732188 

.142738 

.218139 

.652539 

Gama 

Delta 

X  tel 

Hnrtna. 

1.654221 

.359526 

1 

.713518 

1.596867 

.341588 

.704686 

1.530151 

.322133 

■>  m 

,704294 

1.444632 

.301434 

.299734 

,712525 

1.349627 

.286464 

.300420 

,704294 

1.247492 

.262974 

.299351 

.711538 

1.156040 

.243009 

.299734 

702525 

1.072148 

.225374 

.299197 

.701144 

.981559 

.206760 

.299581 

.702(31 

.911912 

.30003? 

.703312 

.039967 

.310039 

.703311 

.781251 

.167502 

.299504 

,701933 

.724876 

.155444 

,299810 

,702721 

Table  II 

(Cont ’ 

d) .  Complete 

-6.800000 

-4.380008 

-2,500000 

0.010000 

2.000000 

4.100000 

6.000000 

3.000000 

10.900000 

12.900000 

14.900000 
16.090090 
18.000000 


Phi 

-6.090000 

-4.000000 

-2.900090 

0.000009 

2.900000 

4.100000 

6.900000 

8.100000 
10.000900 
12.000000 
14.000000 
16.000000 
t8. 000880 


Phi 

-6.000000 

-4.000000 

-2.900000 

0.900000 
2.000000 
4.010000 
6.000000 
8.100000 
10.900000 
12.000011 
14.000000 
16.000000 
18.  HUH 


Phi 

-6.100001 
•4.100100 
-2.000000 
1.000000 
2,011101 
4 . 111888 
6.000091 
B.  900000 
11.000000 
12.000000 
14.000000 
16.100000 
18.0H0I0 


T  Tunnel 
(92.793050 
193.384750 
103.149300 

\m\\ 

103,275741 

103.286160 

103.293120 

103.442460 

103,244510 

(03.414670 

193,432020 

103.494540 


T  Tannel 
103.352230 
104.915460 
104.126570 
104.140470 
104.189090 
114.334951 
(04.345170 
104.328000 
114. 432191 
114.471381 
114.459980 
104.567630 
114.612788 


T  Tannel 
105.776150 
105.960220 
105.970630 
106.213310 
106.241520 
106.335280 
106.380430 
106.488080 
106.581850 
116.585331 
106.621160 
116.734660 
t 04.866621 


T  Tannel 
107.561170 
107.762590 
107.877200 
107.773010 
108.085570 
108.148090 
108.259210 
118.172390 
108.46063! 
108.538081 
108.648160 
108.651630 
108.880841 


T  *ef  jiB* 
104.482S90 
(0^  ,0(5200 
195, '58149 
(85.47063* 
105.149381 
115.973530 
105.286480 
145.268970 
105.341890 
(05.315640 
105.336060 
(05.403140 
105.394390 


T  ref  Dipt 
105.341530 
106.012530 
106.(67020 
(86.164110 
106.172850 
106.324420 
104.295270 
106.382710 
106.414760 

!»:!?$! 

106.615830 

106.674120 


Tctmonbe 
104.533770 
194.T34T28 
'.05,358209 
104,794c64 
(04.822669 
104.816801 
115.108550 
105.(43540 
105. 155210 
t05.076Mi 
105.228150 
1  OS. 191231 
10S. 058960 


T  ref  pipe 
107,336200 
(87.914810 
198.034150 
108.066180 
108.156400 
188.261180 
108.406690 
108.363040 
108.604570 
108.688950 
108.735470 
108.758760 
1H.782BH 


T  ref  pipe 
109.584760 
109.799910 
109,991760 
110.003390 
110.317930 
110.166150 
110,247530 
110.349240 
1(0.378300 
110.570080 
1(1.581700 
110.700820 
110.892560 


Bet  For 
,99S4<8 
09445'; 

.  »06o9c 
J9«*39 
.994806 
.992124 
.994SA8 
.996161 
.994231 
.992652 
.996697 
.993478 
.989686 


TcinCprobt 

(05.677250 

Bet. Far 
.995681 

(05.788060 

.994082 

(05.960050 

.994541 

(06.04(690 

,996743 

(06.062100 

.997974 

(06.263210 

.998382 

(06.164110 

.9945(7 

106.094150 

.992327 

106.190340 

.994042 

lil:» 

m 

106.560470 

.998540 

1M. 338990 

Tcan6pra6e 

107.821621 

Bet. Fac 
.99967? 

107.798360 

.997376 

107.943920 

.997971 

107.874050 

.995677 

107.978840 

.9940(8 

107.958470 

.993153 

108.106900 

.993242 

109.185500 

.995977 

(08.415470 

.9957(3 

108.403790 

.993563 

108.432880 

.993204 

118.496900 

,994171 

1M.S2I19I 

.994134 

Tcintiprabe 
119. 4831 •» 

Bee. Fac 
.998016 

109.683611 

.997741 

110.032460 

1.080791 

119.849320 

.996991 

109.773740 

.995253 

109.924911 

.995278 

1*9.983050 

.994838 

111.087690 

.994877 

111.17778* 

.996083 

110.321181 

.99513S 

111.360870 

.995701 

118.384110 

.993813 

111.6427*1 

.995131 

Complete  Set  of  Reduced  Calibration  Data 
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COEFFICIENTS  FOR  THE  CALIBRATION  SURFACE  STORED  IN  FILE  COEFOX::28 
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TABLE  III.  Pneumatic  Calibration  Coefficients  and  Errors  in 
Surface  Approximation  of  Calibration  Test  Data 


asasaa 


'7FFF*CIENTS  FOR  THE  1  AlIBRATION  SURFACE  STORED  IN  FILE  COEFOR 


\ 

1 

.967377 

-.133312 

-i.249223 

1 

w 

.8294S3 

1,236718 

-1,644683 

7 

w 

,848476 

-1.367321 

7.522S 79 

i 

-.119724 

.257648 

-1.998874 

c 

.184485 

-.8/8578 

-3,3956X7 

5 

-.835893 

,251334 

.381618 

1 

3 

1 

37.A7.Tlt3 

-41,871518 

28.097572 

-8,373269 

12.758778 

-4,679673 

t 


-•91,838888 

-43.778836 

38,816373 

28,396744 

2,119871 

-3,811747 

c 


813.623668 
-339.813868 
-145.381368 
-•9,882 
18.219 
V 286416 
6 


mmCL)  AT  EACH  POINT 
Mach  #  Pitch  #  — 


1 

.847 

-.691 

.481 

.236 

.816 

.139 

.819 

1 


-.849 

.884 

-.233 

.888 

-.232 

-.118 

-.188 


3 

.394 

.395 

.m 

.188 

-.144 

.822 

.235 

3 


4 

5 

6 

n 

9 

7 

18 

M 

12 

f  T 

-.826 

-.718 

-.397 

,S4t 

,747 

-.743 

-.661 

1.657 

-  287 

-.57* 

-.148 

.813 

.981 

-.975 

.185 

-.326 

.125 

,M8 

-'.345 

-.947 

.184 

.336 

.226 

-.271 

-.111 

-.  t44 

.288 

-.383 

.196 

-,1?1 

-.331 

-.446 

-.215 

.931 

.188 

.021 

-.893 

.392 

.178 

-  9ft2 

.115 

,222 

.388 

.195 

-.245 

-.874 

-.253 

.897 

.568 

-.850 

-.078 

■  111 

-.188 

-.187 

.123 

.956 

-.164 

-.164 

.811 

,li» 

-.882 

-.8.16 

.849 

,83? 

-.089 

.851 

-.87 £ 

.888 

-.131 

.  >17 

4 

5 

6 

7 

3 

9 

10 

•1 

12 

'3 

Table  IV.  Temperature  Calibration  Coefficients  and  Errors 
in  Surface  Approximation  of  Calibration  Test  Data 


28 


I 


Pitch  Ancle  Mach  Number 


Table  V.  Probe  Calibration  Verification  Test  Results 
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Figure  7.  Data  Acquisition  System 
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Normalized  Displacement 


APPENDIX  A.  RELATIONSHIP  OF  PROBE  PRESSURE  TO  FLOW 


VELOCITY  AND  ANGLE  (by  R.  P.  Shreeve) 

l-1.  Dimensio.-less  Velocity  and  Mach  Number 

The  velocity  (V)  can  be  made  dimensionless  by  dividing  by 
the  "limiting"  or  "total"  velocity  (Vfc)  which,  for  a  perfect  gas 
is  given  by 


- 


where  Tfc  is  the  total  temperature .  Vfc  is  the  maximum  steady  velo¬ 
city  that  the  flow  can  have  if  the  stagnation  temperature  remains 
unchanged  since,  by  definition 


CpTt  =  CpT  +  V 


and  if  the  flow  is  accelerated  to  the  limit  of  T  0,  the  velocity 
is  that  given  by  Equation  A(l) . 

The  dimensionless  velocity  obtained  by  dividing  the  velocity 
by  its  limiting  value  can  be  viewed  as  the  "Fractional  Velocity" 
or  simply  the  "dimensionless  velocity"  and  is  here  given  the 
symbol  X;  thus 


M3) 


The  relationship  in  Equation  A (2)  can  be  seen  to  give 


A(4) 
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so  that,  using  isentropic  relationships,  the  ratio  of  the  3tatic 
(p)  to  total  pressure  tpfc)  is  given  by 


Pt 


[f]*  -  (l  - 


A(5) 


and  static  Cp)  to  total  density  (pfc)  is  given  by 


A  (6) 


Using  Equation  A (2)  and  the  definition  of  Mach  number  in  a 
perfect  gas,  the  relationship  between  the  dimensionless  velocity 
and  Mach  number  can  be  shown  to  be 


X2 


l^M2 

1  +  M2 


or 


M2 


A  (7) 


A(8) 


Clearly  there  is  a  non-linear  but  one-to-one  relationship 
between  X  and  M,  so  that  if  one  is  known  so,  uniquely,  is  the 
other.  In  compressor  and  turbine  flows  it  is  often  more  desirable 
to  work  with  the  dimensionless  velocity,  X,  since  it  is  directly 
proportional  to  the  velocity  itself  until  the  stagnation  tempera¬ 
ture  is  changed,  as  across  a  rotor.  The  relationships  in  Equa¬ 
tions  A(4),  A(5)  and  A(6)  are  also  simpler  than  the  corresponding 
expressions  involving  Mach  number. 
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In  representing  the  calibration  of  the  probe,  the  dimension¬ 


less  velocity,  rather  than  Mach  number,  was  used. 


A-2.  Probe  Pressure  Coefficients 

The  probe  geometry  is  constant  and  when  used  is  always 
adjusted  so  that  the  probe  is  aligned  (in  yaw)  with  the  flow. 
Therefore,  by  dimensional  analysis,  the  pressure  coefficient, 
Cj^,  defined  as 


CPn 


A  (9) 


where  pn  is  the  pressure  at  sensor  n,  p  and  M  are  the  static 
pressure  and  Mach  number  respectively  of  the  oncoming  flow,  and 
Y  is  the  ratio  of  specific  heats,  can  be  only  a  function  of  the 
Mach  number,  pitch  angle  (<j>)  and  Reynolds  Number  (Rg)  ,  if  a  per¬ 
fect  gas  and  uniform  flow  are  assumed. 

As  a  consequence,  the  difference  between  the  pressure 
coefficients  at  two  sensors,  m  and  n,  defined  as 

CP»m  ’  CPm  '  <=Pn  *  ^2^  i(1° 

must  also  be  only  a  function  of  Re,  M  and  $. 

The  Equation  A (10)  can  be  rewritten  as 


CPmn 


fpm  '  pnl 

fpml 

rpti 

CM 

— 

*o 

9 

# 

la* 
_ / 

e 

p 

i  j 

A  (11) 


where  pt  is  the  stagnation  pressure.  Using  Equation  A (5)  and 
Equation  A(8) ,  Equation  A(ll)  can  be  expressed  as 
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A  (12) 


CPran 


fpm  “  pn^ 

rpm} 

.  Pm 

pt 

1 

•  — 

v 


where  v  is  an  explicit  function  of  the  dimensionless  velocity 
(or  Mach  number) ;  namely, 


1 


Consequently,  using  Equation  A (12),  in  general 


pm  ~  pn 


•m 


Pmn 


A  (14) 


where  the  quantities  on  the  left  hand  side  involve  only  the  sensor 
measurements  and  where  it  is  known  that  Cp^  depends  on  Re,  M  (or 
X)  and  <#> ,  and  v  (X)  is  given  explicitly  by  Equation  A(13)  .  From 
Equation  A (9) ,  Equation  A(5)  and  Equation  A (7)  it  must  follow 
also  that  the  ratio  Pm/Pt  depends  only  on  RQ,  M  and  $ 

In  principle,  from  an  examination  of  Equation  A (14) ,  four 
independent  sensors  are  necessary  in  order  to  calibrate  a  pres¬ 
sure  probe  uniquely  for  variations  in  M,  Rg  and  <P .  The  present 
probe,  when  rotated  to  balance  p2  and  p3  (approximately)  has  only 
three  independent  measurements;  namely,  p^,  p23  and  p^,  where 
P23  is  the  pneumatic  average  of  the  nearly-equal  pressures  p2 
and  P3.  Only  two  independent  parameters  of  the  type  in  Equation 
A(14)  can  therefore  be  written.  They  are  defined  as 


8 


P1  ~  P23 
pl 


and 


5 


pl  ~  p4 
pl 


A(15) 


A  (16) 
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where  an  alternate  for  either  one  is  obtained  in  the  ratio 


5  P1  ~  p< 
*  “  pl  ‘  p; 


A  (17 ) 


However,  the  geometry  of  the  sensors  is  such  that  the 
effect  of  Reynolds  number  on  the  pressure  coefficients  should  be 
extremely  small  in  the  range  of  Mach  number  of  interest,  and  this 
has  been  verified  experimentally  (Reference  6  ) .  Therefore,  in 
principle,  calibration  is  required  to  establish  the  relationships 


S  -  B(X,4) 


A  (18) 


«(X,*) 


A  (19) 


Y  *  y  (X,<j>)  . 


A  (20) 


The  selection  of  B,  S  or  y  is  not  arbitrary  because  it 
involves  a  selection  between  different  sensitivities.  For  example, 
from  Equation  A (14) 

„  _  P1  '  p23  _  Cp12  .  _ 

8 - Pi - u(x) 

Since  CPl2  will  probably  not  be  strongly  sensitive  to  Mach  number, 
and  (Pj/Pt)  will  be  very  close  to  unity  for  moderate  pitch  angles, 

8  is  a  direct  measure  of  the  function  v(X)  and  therefore  of  the 
Mach  number,  with  a  second  order  dependence  on  pitch  angle.  In 
contrast 


pi :  p4 

pl  ’  p23 


48 


will  not  depend  much  on  Mach  number  but,  through  Cpi4,  because  of 
the  orientation  of  the  sensors,  will  be  highly  sensitive  to  pitch 
angle.  The  quantity  y,  given  by 


Y 


P1  -  P4 
P1 


CP14 
(Pl'/Pt)  * 


y 


will  again  be  directly  a  measure  of  Mach  number  but,  because 
of  CP14,  will  also  be  strongly  dependent  on  pitch  angle. 
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APPENDIX  B.  DATA  ACQUISITION  AND  REDUCTION 

B-l.  Data  Acquisition  Using  BASIC  Program  AQCPRB 

The  data  acquisition  was  carried  out  using  a  Hewlett 
Packard  Model  9830A  Calculator  as  controller.  BASIC  language 
program  AQCPRB  was  written  and  stored  on  mass  storage  platter 
PL-001.  A  program  listing  is  given  in  Section  B-3.1. 

The  connections  of  the  instrumentation  shown  in  Figure  6 
to  the  data  acquisition  system  shown  in  Figure  7  are  given  in 
Table  B-I.  The  data  files  containing  the  raw  data  stored  on 
PL-001  are  given  in  Table  B-II. 

B-2 ,  Data  Reduction  Using  FORTRAN  Programs  REST2  and  COEFS 
The  data  reduction  was  carried  out  using  a  Hewlett 
Packard  Model  HP21MX  computer  (Figure  7) .  First,  the  data 
files  given  in  Table  B-II  were  transferred  to  the  HP21MX  using 
the  system  utility  FORTRAN  program  X9830  (stored  on  cartridge 
28)  on  the  HP21MX,  and  BASIC  program  9830X  on  the  HP9830A.  The 
data  file  names  were  retained  on  the  HP21MX. 

Second,  Program  REST2  (cartridge  28)  was  used  to  collect 
the  raw  data  files,  reduce  the  data  to  pressure  coefficients, 
temperatures  and  recovery  factors,  and  then  to  store  the  re¬ 
duced  data  in  a  single  large  data  file  CALNEW  (on  cartridge  28) . 
The  complete  reduced  data  in  CALNEW  are  shown  in  Table  B-III 
and  the  listing  of  program  REST2  is  given  in  Section  B-3.2. 
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The  process  of  obtaining  the  coefficients  of  surface 
approximations  to  the  calibration  data  was  carried  out  using 
program  COEFS  (on  cartridge  28)  which  is  listed  in  Section 
B-3.3.  The  program  reads  the  reduced  data  file  CALNEW  into 
an  array.  The  surface  approximation  procedure  is  carried  out 
for  the  dimensionless  velocity  X,  the  pitch  angle,  or  the 
temperature  recovery  factor  as  required;  one  at  a  time.  Which¬ 
ever  is  required  must  be  specified  in  all  the  statements  which 
carry  an  asterisk  in  column  #73.  (In  the  listing  in  Section 
B.3.3  it  is  the  pitch  angle). 

To  output  the  plot  correctly,  the  statements  with  aster¬ 
isks  in  columns  #73  and  #74  must  be  edited  to  adjust  the  range 
of  the  Z-axis. 

The  order  of  the  approximation  is  input  by  the  operator 
and  can  be  varied  interactively  up  to  an  order  of  6  for  both 
independent  variables.  The  coefficients  for  the  approximations 
and  the  errors  for  each  data  point  (see  Equations  4,  5  and  10) 
are  printed  and  plotted  and  the  operator  is  asked  for  a  name 
for  the  file  in  which  to  store  the  coefficients. 

B-3.  Program  Listings 

The  following  pages  contain; 

B-3.1  Basic  Program  ACQPRB 
B-3. 2  Fortran  Program  REST2 

B-3.3  Fortran  Program  COEFS. 
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10  R£H******AQCPRB**  *********************************** 
20  REM******ON-L INE  DftTft  ACQUISITION  PROGRAM*********** 
30  REM******FOR  ft  COMBINATION  PROBE******************* 

40  REM 

50  REM—  NOTE  :  THIS  PROGRAM  ONLY  ACQUIRES  THE  RAW  BATA 


B-3.1.  Basic  Program  AQCPRB 
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PHI  9ECF 


MACH  Ort&F.R  =  5 


5AMMA  DELTA  X'JEL 
.'tfCH  NUMBER  =  i 


1.52)552 
1. 160879 
1.397381 


4812882 
.0773556 
. 5736720 


1,313777  .3695280 
1.254386  .8658292  , 
i. 170305  .3617003 
1.100437  .0580034  , 
1.325362  . 0S47432 
.9521739  .0504077 
.3744588  . 3465162 
.3034189  .8432724  , 
.7531915  .0407311 
.7008547  .0377656 

hack  number  =  2 

1.535627  .1371644 
1.437562  .1311814 
1.420398  .1252310 
1.359802  .1201603 
1.277915  .1128996 


,1332521 
.1327658 
1320334 
.1327658 
132.' 


132522' 

1325223 

1325223 

1322778 

1327658 

13300  9 

1332521 

1327658 


-6.30000 

-4.30000 

-2.30000 

0.300000 

2.300000 
4.900000 

6.300000 

3.300000 

10.30000 

12.30000 

14.30000 
16.30900 
18.90000 


1.182716 

1.103961 


1050076 

.3978161 


1.012255  .3905588 
.9511302  .1852776 
.8725491  ,0781117 
,3097562  .0728617 
.7628361  .0684725 
7094432  .0643027 
MACH  NUMBER  *  3 

1.560976  .1880941 
t. 507067  .1791348 
1.448028  .1696561 
1.375000  .1616433 
1.289521  .1523746 
1.197158  .1414903 
1.116608  .1326734 
1.329982  .1226227 
.9543859  .1142001 
.3778360  .1056373 
.8124999  .0982455 
.7604167  .0919478 
.7137872  .0859320 
HACK  NUMBER  *  4 


.1763513 
.1768670 
.1768670 
.1761789 
.  176S233 
.1765233 
.1766954 
.1768670 
.1765233 
, 1766954 
.1760063 
.1770385 
1760063 


-6.00000 

-4.00000 

-2.0(000 

0.000000 

2.000000 

4.000000 

6.000000 

8.000000 

10.00000 

12.00000 

14.00000 

16.00000 

18.00000 


.9917194 

9886776 

.9915637 

.9862901 

.9916281 

.9937222 

.9789860 

.9800366 

1.002970 

,9840231 

.9813301 


.8089529 

.8213995 

.7210634 

.7059393 

.7402264 

.6759887 

.6942834 

.7305281 

,6986769 

.7744532 

.7088474 

.7349200 

,6128314 


.2077222  -6.80000 
.2085493  -4.00000 
2078605  -2.00000 


.20' 


0.000000 


2078605  2.000000 
,2079984  4.880008 

-.3HMBH 

.2081361  10.00000 
.2079984  12.00000 
.2084118  14.00000 
.2178615  16.00000 
2182741  18.00000 


98S3598 
.9934516 
9965249 
,991*246 
,99 


1.576882 

1.513386 

1.455272 

1.373206 

1.285714 

1.196172 

1.113924 

1.031496 

.9527559 

.8775511 

.8159127 

.75937S1 

.7138366 


2075264 

1981208 

1878903 

1776145 

,1670249 

,1547484 

,1451374 

,1350634 

.1248147 

,1153392 

,1078667 

.1002563 

.0937518 


2191589 

2190310 

2186459 

2187743 

2191589 

,2192869 

2194147 

,2192869 

,2189025 

,2195426 

2191S89 

,2197979 

,2187743 


-6.00000 

-4.00600 

-2.11000 

0.000000 

2.000000 

4.000*80 

6.000000 

8.900*00 

10.00000 

12.00000 

14.00000 

16.00*00 

18.90000 


.9912542 

.9929143 

.9884121 

.992.3444 

.9945201 

.9861148 

.9880694 


.9968155 

.9908348 

.9955633 

.9943811 

.9965537 

,9952673 

.9903313 

.9918269 

.9911489 

.9939903 

.9882616 

.9921843 

,9845446 


1.595331 

,538058 

.,473545 


.2448943 

.2333929 

.2216220 


t, 394422  ,2088067 
1,300662  .1954783 
1,208719  ,1821775 
1,1)9423  ,1698332 


1.043081 

.9595301 


.1588919 

,1462810 


,6844156  ,1355339 
,3223087  ,1262552 


,7673338 

.7148439 


,1178678 

.1094591 


,2406194 

,24*2852 

.2403966 

,2406194 

.240S079 

.2401736 

.2407386 

.2403966 

2402852 

,2399504 

2403966 

,2402858 

2398386 


1ACH  NUMBER  *  6 

1,609131  ,2783892  .2584729 
1,553336  .26656*3  ,2588692 
1,484055  ,2506939  ,2587782 
1.411968  .2361171  ,2576768 
1.3)5428  .2215773  ,2585720 
1.217540  ,2057917  ,2584729 
1,126554  .1918572  .2578761 
1.048533  .1787716  .2576768 
,9672318  ,1648825  .2578761 
.8901345  .1528812  .2585718 
,8234637  .1419061  .? 587702 
,7731843  .1332672  .2586711 
,7170023  ,1234930  .2579757 
MACH  NUMBER  =  7 

1.64)706  ,3195823 
1.583654  .3039544 


i. 506256 
1.427875 


2881830 

.2704661 


1.329150  .2540321 
1.232625  .2356273 
1.146035  .2185319 
1,058027  .2019354 
,9731286  .1869960 
,8975096  , 1728S99 
,8358778  .1616661 


,7793696 

,7320175 


.1507600 

1427365 


MACH  NUMBER  =  A 


.2796194 

.2795326 

,2797930 

,2797564 

2798797 

,2789234 

2793587 

,2786617 

2785742 

.2790977 

2797930 

,2786617 

2801394 


-6,00000 

-4.00000 

-2,00000 

0.050000 

2.000000 

4.000000 

6.000000 

8.000000 

19.00*00 

12.0**00 

14.00000 

16.00000 

13.00000 


-b. 00000 
-4.00000 
-2.50000 
0.000000 
2.009099 
4.000009 
6.000000 
8.000000 
10.90000 
12.00000 
14.00600 
16.00000 
18.00000 


-6.00000 

-4.00000 

-2.00000 

0.000000 

2.000000 

4.000000 

6.0000(0 

8.000000 

10.00000 

12.00000 

14,00009 

16.00009 

18.900(0 


1.654221 

1,596867 

,530151 

..444632 

1,349627 

1.247492 

1.156(40 

1.072148 

,98)5592 

,9019120 

,8399669 

.7812550 

,7248764 


3595260 

3415878 

,3221327 

3014336 

2864642 

2629735 

2430088 

2253740 

2067599 

1913724 

1787676 

1675020 

1554063 


,3001155 

.3005725 

.3004204 

.2997335 

,3004204 

.2993505 

,2997335 

.2991971 

,2995806 

,3(00392 

301(391 

,2995039 

,2998099 


-6.909(0 

-4.90*00 

-2.1(000 

9.900001 

2.900(00 

4.000*00 

6.900000 

8.000*00 

10.900(9 

12,00*00 

14.00011 

16.00*00 

18.90000 


,9954481 

,  9944553 

.9966950 

.9914392 

,9999056 

.9921237 

.9945685 

.996) 60S 

,9942899 

.9926516 

,9966971 

,9934781 

.9896863 


,9956807 

.9940828 

,9945407 

.9967434 

.9970744 

.9983823 

.9965171 

.9923269 

.9940418 

.9919972 

.9947754 

.9985398 

.9911137 


.9996718 

.9973757 

.9979711 

.9956771 

,996(102 

.9931526 

.9932418 

.9959774 

.9957135 

.9935628 

.9932*35 

.994*710 

.9941336 


,998*165 

.9977407 

1.000791 

,9969914 

.9952534 

.9952784 

.9948375 

.9948773 

.9960829 

.9951346 

.995701* 

.9938134 

,99513*8 


Table  B-ill.  Reduced  Calibration  Data 
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APPENDIX  C.  CHROMEL-CONSTANTAN  THERMOCOUPLE 
OUTPUT  APPROXIMATION 

The  relationship  between  the  millivolt  output  of  the 
Type  E  (Chromel-Constantan)  thermocouple  element  and  the  tem¬ 
perature  is  given  by  the  manufacturer  (Omega  Engineering, 
Inc.,  P.0.  Box  4047,  Stamford,  Connecticut,  06907)  from  NBS 
Standard  tables.  For  a  range  of  0°  to  170°  Fahrenheit  the 
temperature  as  a  function  of  the  voltage  was  approximated  by 
a  second  order  polynomial  using  the  method  of  least  squares. 
Figure  Cl  shows  the  approximation  together  with  the  coeffi¬ 
cients  for  the  curve.  The  second  order  approximation  in 
Figure  Cl  was  used  in  the  data  reduction  programs  for  the 
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Figure  Cl.  Chromel-Constantan  Thermocouple 
Voltage  Output  vs  Temperature 


APPENDIX  D.  PNEUMATIC  AVERAGING 
OF  UNSTEADY  PRESSURES 


Other  investigators  have  questioned  whether  a  pneumatic 
probe  measures  the  correct  time  average  value  of  the  impact 
pressure  in  an  unsteady  flow  (References  7  and  8) .  It  is  clear 
that  the  possible  error  in  the  probe  measurement  will  depend  on 
the  relationship  between  the  average  pressure  level  and  the 
unsteady  departures  therefrom.  Weyer,  in  Reference  7,  has 
measured  the  magnitude  of  the  error  for  a  particular  wave  form. 
In  order  to  get  an  idea  of  the  error  which  might  be  involved  if 
the  measurement  of  PI  was  taken  to  be  the  time-averaged  total 
pressure,  the  unsteady  pressure  distribution  in  the  measuring 
plane  of  the  probe  was  observed  in  the  compressor.  Figure  Dl 
shows  the  pressure  distribution  measured  behind  the  compressor 
rotor  using  a  simple  impact  pressure  probe  incorporating  a 
Kulite  semiconductor  transducer  at  the  tip.  Rotor  rotational 
speed  was  15,200  rpm,  resulting  in  a  (blade-passing)  frequency 
of  4.56  Khz  in  the  pressure  variation.  Using  the  notation  of 
Reference  7,  the  characteristics  of  the  signal  were  determined 
to  be  the  following: 


Time  average  value  of  the  pressure,  P1  *  434"  H20 
Peak-to-peak  amplitude  of  the  pressures. 


(p««  -  w  ■  2P1« 


16"  HjO 


Ratio  of  wave  amplitude  to  mean  value. 


0.037 
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Figure  D2  (from  Reference  7)  3hows  the  relationship  between 
the  error  which  was  measured  for  a  particular  pneumatic  probe, 
the  frequency,  the  pressure  level  and  the  amplitude  of  the  pres¬ 
sure  changes  for  the  particular  wave  shape  indicated  in  the 
figure.  It  is  observed  that  the  frequency  involved  in  the 
present  measurements  exceeds  the  range  given  in.  Figure  D2,  and 
that  the  value  of  221^/?^  was  also  very  much  smaller. 

Similar  comparisons  between  the  pressure  readings  P2  and 
P3  and  corresponding  Kulite  probe  measurements  showed  that  for 
these  sensors  2Pla/P^  was  even  smaller,  thus  any  error  in  the 
averaging  of  P2  and  P3  should  be  negligible. 

An  equivalent  Kulite  probe  for  P4  did  not  exist,  so  that 
a  comparison  for  this  sensor  could  not  be  made.  However,  since 
the  tube  for  P4  is  inclined  at  an  angle  to  the  flow  comparable 
to  those  of  P2  and  P3,  it  can  reasonably  be  assumed  that  the 
conclusions  are  similar. 

For  operating  conditions  at  higher  speeds  and  flow  rates, 
the  above  conclusions  must  be  reexamined  for  the  time-varying 
rotor  exit  conditions  then  measured.  Whereas  the  data  in  Figure 
D2  suggest  that  the  error  becomes  small  at  higher  frequencies, 
it  was  shown  in  Reference  9  that  measurable  errors  could  occur 
at  3  -  8  Khz  in  particular  pneumatic  systems  if  the  wave  form 
was  sufficiently  extreme  (implying  large  values  of  the  parameter 

2PVF1>  • 
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